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Abstract 

cn" 

A discussion is given of the quantisation of a physical system with finite degrees of freedom subject 
to a Hamiltonian constraint by treating time as a constrained classical variable interacting with an 
unconstrained quantum state. This leads to a quantisation scheme that yields a Schrodinger-type 
equation which is in general nonlinear in evolution. Nevertheless it is compatible with a probabilistic 
interpretation of quantum mechanics and in particular the construction of a Hilbert space with a 
Euclidean norm is possible. The new scheme is applied to the quantisation of a Friedmann universe 
with a massive scalar field whose dynamical behaviour is investigated numerically. 
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^ ■ 1 Introduction 
o 

*q I An elegant way of formulating classical particle dynamics is via an action principle that is invariant 
under a general change of time parameterisation. Such a formulation is called a parameterised theory. 
In non-relativistic particle dynamics a parametric description is obtained by regarding Newtonian time 
of the particle as a dynamical variable depending on a newly introduced "parameter time" pQ. As 
this choice is arbitrary the corresponding Lagrangian admits a symmetry under re-parameterisation 
^ . of time, leading to a constraining relation called the Hamiltonian constraint. Consequently not all 
canonical variables (including the original Newtonian time and its conjugate momentum) can be given 
arbitrary initial data. Rather, the Hamiltonian constraint must be satisfied initially which will be 
preserved under the canonical equations of motion. This is to be expected since parameterising time 
cannot change physical kinematics, the apparent extra degree of freedom associated with time must 
be subject to a constraint. 

In relativity the parameterised theory is ideally suited for formulating particle dynamics where 
time is treated on an equal footing with other coordinates. This is especially the case for a particle 
moving in a (pseudo) Riemannian geometry The canonical description of a parameterised rela- 
tivistic particle is similar in structure to that of a non-relativistic particle. However the constrained 
Hamiltonian is now quadratic in the canonical momenta of both position and time coordinates. This 
is due to the (general) covariance of relativistic kinematics and can be analysed within the established 
methodology in classical theory. However this appears to present a number of challenges for the 
quantisation of a relativistic particle compatible with standard interpretations. As is well-known, a 
straightforward application of canonical quantisation results in a Klein-Gordon equation. Unfortu- 
nately its solutions belong to a function space with an indefinite norm and cannot be interpreted as 
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wavefunctions describing a probability amplitude in the usual manner. Although this difficulty can 
be circumvented upon "second quantisation" in the context of relativistic quantum field theory, the- 
oretical interest remains in search for a satisfactory direct quantisation of a relativistic particle. One 
important reason for this is that a functional form of relativistic particle dynamics naturally arises 
from canonical general relativity 3 . Furthermore systems with finite degrees of freedom subject to a 
constrained Hamiltonian formally resembling that of relativistic particle dynamics in curved spacetime 
often feature in various cosmological models with or without coupled matter fields. These models are 
normally constructed by imposing high symmetry on spacetime geometry and matter distribution in 
order to gain insight into aspects of an underlying full theory. Note that in these models the "metric" 
refers to the (reduced) "supermetric" (of spacetime metric) and the "relativistic particle" represents 
a collective gravitational and matter field point. 

In |3| a "square-root Hamiltonian" method was investigated in an attempt to overcome the prob- 
lems mentioned with the "Klein-Gordon approach" . Given a choice of time coordinate it is possible to 
construct a wave equation that is first order in time and possesses a "Schrodinger structure" . While 
in some cases this structure guarantees the existence of a Hilbert space with a Euclidean norm it may 
suffer from other drawbacks including the violation of Hermiticity of the Hamiltonian operator if the 
square-root of a negative eigenvalue must be taken through a spectral analysis [3]. 

Motivated by these circumstances an alternative canonical type of quantisation for a physical 
system with finite degrees of freedom subject to a Hamiltonian constraint is considered by treating 
time as a constrained classical variable interacting with an unconstrained quantum state. Following 
a brief review of the basic parameterised theory of a non-relativistic particle in curved space and its 
Dirac quantisation in section a new "parametric quantisation" is introduced. For non-relativistic 
particle dynamics, this new scheme is shown to reproduce results based on Dirac quantisation. In 
section |3 parametric quantisation of a relativistic particle in curved spacetime is performed for a 
choice of time-slicing. In this case a "nonlinear Schrodinger equation" may be derived that is locally 
equivalent to the full set of evolution equations from parametric quantisation under certain conditions. 
In section 0] the theory is applied in the context of a simple Friedmann cosmological model with the 
resulting equations analysed numerically. Concluding remarks are made in section |S] where some 
prospects for future work are also discussed. 



2 Classical and quantum dynamics of a non-relativistic 
particle 



Consider the motion of a non-relativistic particle subject to a potential in a Riemannian n-space. 
Denote by r the Newtonian time and introduce the coordinates (q a ) (0 < a < n with q° := r). 
The time-dependent spatial metric components are denoted by 7 a fc = jabiQ 01 ) = lab{<f-, T )- Here and 
throughout Greek letters such as a, (3 denote spacetime indices with range 0, 1, • • • n, whereas Latin 
letters such as a, b denote spatial indices with range 1, 2, • • • n. The particle motion is described by 
the trajectory (q a (t)) = (q a (t),T(t)) parameterised by a parameter time t. The choice of t is related 
to a positive function N = N(t) that sets the parameterisation gauge and enters into the generic form 
of the Hamiltonian 



that generates the dynamics of the particle for some 7i. For a non-relativistic particle with canonical 
momenta w := po and p a conjugate to r and q a respectively, H is linear in zu and quadratic in p a as 
given by 
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where V(q a ) = V(q c ,r) is a time-varying potential and the particle's mass has been normalised to 
unity. Based on an action principle using the Lagrangian 

L = q a Pa - H (3) 
where ' = 4?, the equations of motion can be derived in the following canonical form: 

-»%■ 

In addition the Hamiltonian constraint 

n = o (6) 

arises from variations with respect to N . Provided a set of initial conditions for the canonical variables 
q a and p a compatible with © is given, this constraint will be preserved by the evolution of the 
canonical variables satisfying (jlj) and Q for a choice of (positive) N(t). Clearly the particle trajectory 
is invariant under re-parameterisation: t — > t' , N — > N' := g-piV. 

A quantum theory of this particles may be established by applying Dirac quantisation where, in the 
Schrodinger picture, classical constraints are "turned into" operators that annihilate wavefunctions 
describing physical states. Following this procedure, the operator 
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H D {q a ,Pa,T,m) := w + -T^Pal^l^PbT^ + y (7) 

where 7 := det(7 a fe) is constructed by substituting p a — ► p a := —1^^ in Ti and factor ordering so 
that the resulting term quadratic in p a is proportional to the Laplacian operator on a wavefunction 
^(q a , t) of weight \ with respect to the metric 7 a fe. Such a weighting choice is very useful in simplifying 
mathematical expressions when dealing with the wavefunction of a quantum particle moving in a time 
dependent spatial metric. 

On Dirac quantisation the classical constraint © becomes the quantum constraint 

H D ^ = (8) 

on any physical state Although Q does not explicitly involve the parameter time t it is natural 
to interpret r as the evolution parameter for wavefunctions. To see this more clearly we may rewrite 
|H1) in the form 

J^ = 7^ (9) 

OT 

where 

U :=- 7~halh ab Pbl~^ + V (10) 

plays the role of the "true Hamiltonian" as in conventional Schrodinger equations. 

By treating r as the evolution parameter the positive-definite inner product of two wavefunctions 
\&i and ^2 at equal time r is naturally defined as 

(#i,* 2 ) := Jn*2d n q. (11) 

Using this definition the operator 7i is self-adjoint and the standard probabilistic interpretation of 
quantum mechanics applies. In particular if a wavefunction ^ is initially normalised according to 
^) = 1 this normalisation will be maintained under JHJ) and therefore the conservation of the total 
probability holds. 
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The derivation of the non-relativistic Schrodinger equation (jSJ) above is often regarded as support 
for Dirac quantisation. However it is worthy noting that this derivation crucially relies on the fact 
that Ti is linear in the momentum w conjugate to the time variable r. For a relativistic particle, this 
is no longer the case, since the Hamiltonian will be quadratic in all momenta instead. 

It is a purpose of this paper to point out the possibility of deriving @ based on an alternative 
quantisation scheme in which the substitution w —* — l-Mp is not made and the explicit parameter time 
(t) dependence will be retained. The guiding principle is that one of the dynamical variables of a 
parameterised system should play the role of time. In the case of non-relativistic particle dynamics, 
r has been chosen to be this variable. As a time variable it is of "classical nature" but need not be 
treated merely as a parameter. Indeed in classical parameterised theory such a time variable is treated 
as a dynamical variable whose conjugate momentum is subject to a Hamiltonian constraint without 
changing the underlying physical degrees of freedom. In view of this we seek to formulate an analogous 
quantisation procedure where the selected time variable is not quantised but treated as a constrained 
classical variable. The coupling of the classical time variable with other quantised unconstrained 
variables described by a wavefunction will be formulated subject to certain requirements to validate 
essential physical interpretations. These requirements include the existence of a positive-definite norm 
for the Hilbert space of the wavefunctions and the reduction to Q, © and © in the classical limit. 

To achieve this consider the operator 

H(q a ,p a ,T,m) :=zu + H (12) 

where 7i takes the form as given in (|10j) . The operator Ti is constructed in the way similar to that of 
7Y_d with the substitution p a — > p a = — 1-^ in 7i but r, w now remain as classical functions of t. In 

contrast to Dirac quantisation where TLd is used in constructing a quantum constraint equation ©, 
the operator Ti will be used to formulate a Schrodinger-type equation 

i^V = NHiP (13) 

for a wavefunction ip(q a , t) of weight ^ with respect to the metric 7^. Nonetheless, like (|llj) . the inner 
product of two wavefunctions ip\ and ip2 can be defined to be 

(Vi.Va) == / ^2 d n q (14) 



under which the operator TL is self-adjoint. It is readily seen that ()13j) maintains the normalisation of 
the wavefunction (ip,ip) for any r(t),tu(t). This makes it possible to consistently impose the special 
normalisation (tl>,ifj) = 1 and then interpret tp(q a ,t) as the probability amplitude for measurements 
of the particle position (q a ) on a spatial hypersurface at time r{t). For general normalisation, the 
expectation value of an operator O = 0(q a ,p a ,T,zu) may be defined as: 



°> ! -W (15) 



in a standard manner. Sometimes it is useful to denote this expectation value by to explicitly 

indicate that it is evaluated with respect to the wavefunction tft. Using 1)1 3|) and 1)15)1 . the derivative 
of the expectation value ( O ) with respect to the parameter time t can be shown to be 




0,H 



(16) 
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where [ , ] denotes a commutator as usual. The succinctness of this expression benefits from the use 
of weight ^ wavefunctions and their operators. Formulae Q15|) and (|16j) will now be used to construct 
the governing equations for r(t) and w(t). Guided by canonical equations © the equations 

(17) 

are adopted in a "semi-classical" fashion to allow the coupling between variables r, xu and wavefunction 
-0 such that the expected classical limit may be recovered. It can now be observed that, using ((TH|) . 

TC^ is a constant of motion under IjlMj) and (|17|). That is, 

s (*>-"■ (18) 

Therefore the condition 

n) = o (19) 



may be consistently imposed and regarded as the quantum analogue of the classical Hamiltonian 
constraint Q. This condition may be considered as a restriction on the initial classical variables r 
and w for any given arbitrary initial wavefunction if). 

For a non-relativistic particle the quantum description using (|13j) , (|17l) and (|19() is in fact equivalent 
to the single Schrodinger-type equation Q. To see this first note that in this case f = N by using 
1)12(1 and (|17|). Adopt the gauge N = 1 so that we may set t = t. With this choice it is clear at once 
that ifj(q a ,T) is unitarily related to *$>(q a ,T) satisfying Q via 

^(q a : T ) = eX p J OT(r')dr^ ip(q a , r). (20) 

Given a generic classical parameterised theory the procedure leading to equations (|13j) . I|17() 
and (|19|) will be referred to as "parametric quantisation". Although both Dirac quantisation and 
parametric quantisation give rise to equivalent results for non-relativistic particle dynamics, these 
two schemes are in general inequivalent, as will be demonstrated in the next section. 



3 Parametric quantisation of a relativistic particle 

Let us now consider the motion of a particle in an (n + l)-dimensional pseudo-Riemannian manifold 
M. with coordinates (q a ), (0 < a < n), in which the metric has components 7 a ^ = 7 a /3((/ /x ) with 
signature (—,+,+,•••). The classical Hamiltonian for this particle H(N, q a ,p a ) = NTi. with trajec- 
tory coordinates q a (t) and canonical momenta p a (t) parameterised by t associated with a positive 
gauge function N(t) subject to a potential V(q^) is specified by 

n(q a ,p a ) = l^PaPp + v. (21) 

In order to follow closely the methods in the preceding section, the manifold Ai is assumed to admit 
a global time foliation so that (q a ) forms a single Gaussian coordinate chart |7j in which 

Too = -1, 70a = 0. (22) 

With such a time-slicing we denote by r(t) := q°(t) and 7 := det(7 a b) where 1 < a, b < n as in the 
non-relativistic particle case. Thus (|21j) becomes 

H = -^ 2 + ^ 1 ab p aPb + V (23) 
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which generates canonical equations of motion 

^ = N™, ^ = -N™ (24) 
dt dpa dt 8q a { ' 

dr „ T cW dm ^ T dTL 

dt aw dt or 

subject to the Hamiltonian constraint 

H = 0. (26) 
bince U in (|231) resembles we now introduce the operator 

ft := -^w 2 + U (27) 

where 7i takes the form of (j!0|) and is self-adjoint with respect to the inner product 

■= I i>l^d n q (28) 



of any two wavefunctions tpk = 4 , k{Q a , i) (of weight |) with respect to the spatial metric j a b- As in the 
preceding section, parametric quantisation now yields the evolution equations for the wavefunction 
ijj(q a ,t) and constrained classical variables r(t),w(t) as follows: 

i^V = NHif> (29) 



H) = (31) 

where the definition of expectation values follows from (|15jl. As in the non-relativistic case, the nor- 
malisation (ip, ip) = 1 can be consistently imposed and the interpretation of if)(q a , t) as the probability 
amplitude for measurements of the particle position (q a ) on a spatial hypersurface at time r(t) can 
be made. 

It is also possible (at least locally) to reduce (|29|). (|30|) and (|3Tj) for the relativistic case to a single 
wave equation. From (|2*3")) . (|3U|) and (|3*Tj) it follows that f = —Nw where 



m = ±j2(H). (32) 



Hence in regions where > and w < one can explicitly eliminate and set t = r by choosing 
the gauge 

TV = (33) 

w 

and introducing 

*(q a , r) := exp f i Jw^dA ^(q a ,r) (34) 
where ip(q a ,r) satisfies fll?) with (|33J). Since (ft) = (h)^, © becomes 
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We are therefore led to a nonlinear integro-partial differential equation describing the quantum evo- 
lution of a relativistic particle under a chosen time-slicing. Although nonlinear in time, (|35|) does 
not present difficulties in the probabilistic interpretation of its solutions. It follows from the scaling 
invariance of ()35|) . i.e. if ^ is a solution to Q35|) so is C^> for any complex constant C and hence both 
^ and C^> represent the same physical state [HIE]- 



4 Friedmann universe with a massive scalar field 

It is not hard to see that the Gaussian coordinate condition {700 = — 1, 7oa = 0} used in the parametric 
quantisation of a relativistic particle can be extended to the condition {700 = 7oo( r )>70a = 0} by 
changing the coordinate time r — > f T \/— 70o(t0 dr' for some negative 700 (t). With this coordinate 
condition one can repeat the derivations of (|29j) . (|30|) and (|31|) using the operator 

n = ^7 ow 2 + n (36) 

with the same forms of Ti and inner product ( , ) as defined in ()10j) and respectively. 

This slightly generalised coordinate condition makes it easier for the application of the parametric 
quantisation to a simple cosmological model for n = 1, where a Friedmann universe filled with a 
massive scalar field is analysed. 1 In terms of the lapse function N(t) and scale factor R(t) as functions 
of the time coordinate t, the Robertson- Walker metric reads 

g = -N 2 dt 2 + R 2 a (37) 

where a is the standard metric on the homogeneous and isotropic 3-space of constant curvature K, 
with K = 1,0,-1 corresponding to the closed, flat and open cases respectively. The dynamics of 
this geometry coupled to a uniformly distributed but time-dependent scalar field (f)(t) of mass m is 
described by the Lagrangian [I]: 

L = --R 2 + —ft - —R 3 <f> 2 + QNKR (38) 

where ' = 4z. The corresponding action f Ldt is manifestly invariant under time re-parameterisation: 
t -> t', N -» N' := §rN. The scale factor R will now be chosen as the (intrinsic) time variable. 
Accordingly we denote by q° = R^q 1 = (ft with the corresponding conjugate momenta given by 

dL 12R. „ 
Po = ^ = — zrR =■ n (39) 
8R N y 1 

9L R 3 ; 

Pi = ^ = W 4>=: P (40) 

where symbols II and p have been introduced for conciseness in the subsequent discussions. It follows 
that the Hamiltonian H = q a p a — L = NTC where a = 0, 1 with Ti having the form of (|21|) in terms 
of the nonzero metric components and potential as follows: 

700 = -12R, 711 = R 3 (41) 
2 

V = ^-R 3 (i) 2 - 6KR. (42) 
By substituting p — > p := — igj into H we obtain the operator 



1 Units of c = H = 16nG = 1 are adopted. 
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which is to be applied to a wavefunction ip((f), t) (of weight ^ with respect to the 1-dimensional metric 
711 = R 3 ) whose norm at any time t is given by 



For simplicity the normalisation condition 



|^(<M)| 2 d0. 



(44) 



(45) 



will be imposed on the wavefunction tp. On parametric quantisation the evolution equations for the 
wavefunction ip and constrained classical variables R(t) and n(i) are then formulated in accordance 
with ijSnjl and ifTTfl as follows: 

i-^ = iV7^V (46) 



1 dR 
N~dt 

i dn 

iV~d7 



dH F 

dn 



n 

12E 



OR 



3 a 2 
2« 4 a</> 



+ 



3m 



2 j.2 



n 2 



2AR 2 



+ 6K 



H F ) = 0. 



(47) 
(48) 
(49) 



To proceed a (time-dependent) functional basis for ip will be chosen in order to reduce the above 
integro-partial differential equations into a system of nonlinear ordinary differential equations of in- 
finite dimensions. This system may be truncated to facilitate numerical simulation. To this end 

consider 

M0,R) ■= {mR Z ) l, Hk{^&4>) (50) 

for integers k > where 

&(x) :=2- fc / 2 vr- 1 /4 (fc!) -i/2 e - 2 /2 Ffc(2;) 
in terms of the Hermite polynomials H~k(x). The functions ip^ satisfy the eigen equations 



i a 2 m 2 ,A , / 1. 
2R i d6 2 2 Y J V k V 2 ' 



and orthonormality 

/oo 
i(>j(<j>,R)ipk(<j>,R)<i(f> = $jk 
-oo 

for any positive i?. They are used to expand the normalised wavefunction: 

oo 

i>(<f>,t) = J2 c k(t)M<f>,R(t)) 

fc=0 



(51) 
(52) 
(53) 

(54) 



with complex coefficients c^it) satisfying XlfcLo \ c k\ 2 = 1- The time-dependence of these coefficients 



is to be determined by dynamics. Using the expansion (|54|) it follows that 



dip 



ICfc + 



livn 

16i? 2 



(^(fc-l)ct-2 - V(fc + 2)(fc + l) Cfe+2 ) 



and 



fc=0 



m il 



n 2 

2AR 



6KR 



(55) 
(56) 
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With the aide of ([52(1 . (|53j) . (|55|) and (j56|) . we obtain the following coupled differential equations for 
Cjfc(i) (fc > 0), R(t) and II(i) by substituting (EU into (HH|) -(H!H): 



Cfc 

iv =1 


r -n-9 

n 2 


24i? 


J? 


n 


iV ~ ~ 


12R 


ri 

N ~ ~ 


oo r 

E- 



+ 6KR 



k -\ — | m 
2 



n 



16i? 2 L 



V(A; + 2)(fc + l)c fe+2 - \A( fc - 1) 



n 2 



fc=0 



R 



24R 2 



+ 6K 



(57) 
(58) 
(59) 



together with the constraint 



fc=o ^ ' 



(60) 



derived from (|49|). and the normalisation condition 



k=0 



1. 



(61) 



Based on these equations numerical simulations are performed for different choices of the curvature 
parameter K to explore the dynamical behaviour of the model. With the mass of the scalar m = 10 
and gauge N = 1 the initial scale factor R(t) at t = is set to 0.5. The initial wave profile 0) 
is chosen to be gaussian by setting co(0) = 1 with Cfc(0) = for k > 1, corresponding to an initially 
"ground state" -like wavefunction compatible with the normalisation condition (|61j) . The initial value 
for the momentum n(0) is then determined by solving the constraint (|49|) at t = with Cfc(0),i?(0) 
substituted with their initial values above. This yields two signs for 11(0) equivalent to integrations 
forward or backward in time t (with a negative or positive n(0) respectively). With these initial data 
and parameters, a truncated k index range < k < 60 is found to produce numerical results with 
satisfactory accuracy. 

For an open Friedmann universe with K = —1 the simulated R(t), H(t) and |^>(0, t)\ 2 are displayed 
in figures for —0.15 < t < 3.0. The wavefunction ip is evaluated using the expression (|50|) in terms 
of the numerical solutions for c^(t). Figure ^ clearly exhibits an expanding universe with the scale 
factor R growing in t and approaching a singularity (R — > 0) with decreasing t. Since R does not 
change the sign of H(t) remains negative as show in figure [21 Figure shows that the probability 
density \ip\ 2 with respect to the measure d(p takes the shape of a single packet. As the universe 
evolves this packet becomes narrower and its height seems to reduce to zero as t decreases to where 

2 



R — ► 0. Despite this, the alternative probability density 



R 



3 

~4ip 



(corresponding to the use of the 



weight wavefunction R it/) with respect to the measure Rzdcft) exhibits a more steady behaviour 

3 2 



as a function of Rt(f) remains the shape 



in its evolution. As shown in figure 0] the profile of R 
of a single packet whose width simply oscillates by either increasing t or decreasing t even to where 
R — > 0. The regular behaviour of R~±ip suggests that the limiting value of this function may serve 
to define the initial quantum state of the Friedmann universe as R — > 0. 

Similar behaviours for R, U and ip can be observed from simulation results for a flat Friedmann 
universe with K = 0. The numerical solutions for R(t), n(t) and \ip(<f>, t)\ 2 for a closed Friedmann 
universe with K = 1 are displayed in figures [5HZI for —0.23 < t < 2.27. Singularities are encountered 
as R — > at both bounds of the numerical integration as shown in figure El Notably in figure n 
changes sign when the local maximum of R takes place. Nevertheless this does not cause any problem 
for the evolution of the system. As in the K = — 1 case the probability density \ip\ 2 with respect to the 
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measure d(p has the shape of a single packet. As t approaches the bounds of integration where R 

3 2 q 







the height of this packet also tends to zero (figure EJ) whereas 



R~il> 



the shape of a single packet with a steadily oscillating width throug 



as a function of R 2 remains 
lout the integration domain as 



shown in figure|Sl This provides further support for the use the limiting value of R ^ as a function 
of R2(p to identify the initial (or final) quantum state of the Friedmann universe as R — > 0. 



5 Concluding remarks 

A quantisation scheme for a physical system subject to a Hamiltonian constraint is presented based 
on canonical quantisation, nonlinear quantum mechanics, couplings between classical and quantum 
variables and analysis of the true degrees of freedom in geometrodynamics 0. Motivated by the 
need for a clarified role played by "time" in the ongoing efforts to ultimately quantise gravity, the 
idea of treating time as a constrained classical variable is put forward in this paper. Instead of 
regarding time as an evolution parameter, the proposed idea treats it as a dynamical variable while 
maintaining its classical nature. It does not follow that the physical degrees of freedom will increase 
by giving time a dynamical status since the time variable and its canonical momentum are subject 
to the (quantised) Hamiltonian constraint. This does imply that if the Hamiltonian constraint only 
restricts the choice of classical time and its canonical momentum, then the quantum degrees of freedom 
can be regarded as unconstrained, thereby enabling much of the standard quantisation technique to 
continue to apply. Of course the issue of the preferred choice of time remains just as in the square-root 
Hamiltonian approach. However it is worth emphasising that the proposed parametric quantisation 
method is free from problems associated with taking square roots of negative eigenvalues arising 
from a spectral analysis. Therefore the choice of a time coordinate is less restrictive for parametric 
quantisation. Whether there exists a naturally preferred time variable given a classical parameterised 
theory or possibility of establishing equivalence for different choices of time within the framework of 
parametric quantisation stays unaddressed in this paper and will be considered in a later publication. 
Finally in extending the present work to the quantisation of gravity it is conceivable that one of the 
traditional functional degrees of freedom in geometrodynamics should be treated as a constrained 
classical variable within the framework set up by this paper. York's "extrinsic time" ^U] is amongst 
potential candidates for such a treatment. In this respect it is of interest to compare with the approach 
by The progress of these further researches will be reported elsewhere. 
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Figure 1: Simulated R(t) for an open Friedmann universe (K = — 1). An expanding universe with the 
scale factor R growing in t is clearly exhibited. This function approaches zero by decreasing t where the 
Robert son- Walker spacetime becomes singular. 
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Figure 3: Behaviour of \ip((f>,t)\ based on simulated 4>(4>,t) for an open Friedmann universe (K = — 1). 
This expression as a function of at a constant t takes the shape of a single packet. As the universe 
evolves this packet becomes narrower and its height reduces to zero as t decreases to where R — > 0. 




Figure 4: Behaviour of R ±iJj based on simulated ip(<f>,t) and R(t) for an open Friedmann universe 

(K = —1). The profile of this expression as a function of R^cj) remains the shape of a single packet whose 
width simply oscillates by either increasing t or decreasing t even to where R — > 0. 
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Figure 5: Simulated R(t) for a closed Friedmann universe (K — 1). Scenarios of Big Bang and Big Crunch 
as R — > are exhibited as in classical cases. 
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Figure 6: Simulated n (t) for a closed Friedmann universe (K = 1) This function changes sign when the 
local maximum of R takes place. This does not affect the continuous evolution of the system. 
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Figure 7: Behaviour of \ip((f>,t)\ based on simulated ip(<f>,t) for a closed Friedmann universe (K = 1). 
This expression as a function of (ft at a constant t takes the shape of a single packet. As t approaches the 
bounds of integration where R — > the height of this packet also tends to zero. 
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Figure 8: Behaviour of R ±ip based on simulated ip(<j),t) and R(t) for a closed Friedmann universe 

(K = 1). The profile of this expression as a function of R^4> remains the shape of a single packet whose 
width simply oscillates by either increasing t or decreasing t even to where R — > 0. 
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